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Abstract: Directional detection of non-baryonic Dark Matter requires 3D reconstruction of low 
energy nuclear recoils tracks. A gaseous micro-TPC matrix, filled with either ^He, CF4 or C4H10 
has been developed within the MIMAC project. A dedicated acquisition electronics and a real time 
track reconstruction software have been developed to monitor a 512 channel prototype. This auto- 
triggered electronic uses embedded processing to reduce the data transfer to its useful part only, 

1. e. decoded coordinates of hit tracks and corresponding energy measurements. An acquisition 
software with on-line monitoring and 3D track reconstruction is also presented. 

Keywords: Electronic detector readout concepts (gas, liquid); Particle tracking detectors 
(Gaseous detectors). 
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1. Introduction 

Directional detection of dark matter is known to be a promising search strategy of galactic Dark 
Matter [U, 0]). Recent studies have shown that, within the framework of dedicated statistical data 
analysis, a low exposure directional detector could lead either to a high significance discovery of 
galactic Dark Matter [0, 3] or to a conclusive exclusion [Q]. 

A gaseous micro-TPC matrix, filled with either ^He, CF4 or C4H10 has been developed within 
the Micro TPC MAtrix of Chambers (MIMAC) project [g]- To demonstrate the relevance of the 
concept, specific front-end ASIC and a dedicated acquisition electronic were developed in order 
to equip a prototype detector featuring an anode of 10.85 x 10.85 cm^ where 2 x 256 strips are 
monitored. This auto-triggered acquisition electronic uses embedded processing to reduce data 
transfer to its useful part only, i.e. decoded coordinates of hit tracks and corresponding energy 
measurements. To be fully exploited, an acquisition software with on-line monitoring and track 
reconstruction has been written. 

2. MIMAC detector readout principle 

As shown in figure [I], the MIMAC prototype pTFC is composed of a pixeUzed anode featuring 2 
orthogonal series of 256 strips of pixels (X and Y) [0] and a micromesh grid defining the delimita- 
tion between the amplification (grid to anode) and the drift space (cathode to grid). Each strip of 
pixels is monitored by a current preamplifier and the fired pixel coordinate is obtained by using the 
coincidence between the X and Y strips (the pixel pitch is 424 \im). A coincidence is defined as 
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Figure 1. Schematic of the MIMAC micro-TPC using a micromegas composed of a pixeHzed anode featur- 
ing 2 orthogonal series of 256 strips of pixels and a micromesh grid defining the delimitation between the 
amplification (grid to anode) and the drift space (cathode to grid). 



having at least one strip of pixels fired in each direction (X, Y) at the same sampling time. The ion- 
ization energy of the recoil energy is obtained by instrumenting the micromesh grid with a Charge 
Sensitive Preamplifier (CSP). 

As illustrated in figure 0, the coordinates in the anode plane (X, Y) are reconstructed by col- 
lecting primary electrons produced in the drift region. Knowing the electron drift velocity, the 
third dimension (Z) is obtained by samphng the anode signal every 20 ns. Note, that due to the 
multiplexed readout of the anode, each time slice picture is rectangular. 




Figure 2. The coordinates in the anode plane (X, Y) are reconstructed by collecting primary electrons 
produced in the drift region. Knowing the electron drift velocity, the third dimension (Z) is obtained by 
sampling the anode signal every 20 ns. Due to the multiplexed readout of the anode, each time slice picture 
is rectangular. 
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3. Front end ASIC 
3.1 Requirements 

In the early stage of the project it was decided to design an ASIC in order to be able to fulfill the 
final objective, which is to equip about 2500 chambers of 1024 strips of pixels (512+512). This 
minimizes space requirement and power demand while allowing cost reduction on a large scale. 
After going through a first prototype phase, 16 channels ASICs equipping a 2 x 96 strips of pixels 
chamber |M]> a 64 channel version was designed [Q]. This was determined to be a good balance 
between integration scale on one side and complexity, fabrication yield and available packages on 
the other side. 

To be able to recover the third coordinate (Z) of the track, a fast switching current compai^ator 
having a threshold as low as 200 nA must be designed in order to have a precise time over threshold 
measurement of each current preamplifier output. This requirement is driven by the worst case 
where the recoil energy is as low as 500 eV in a chamber having its gain limited to 3000 and 
where the diffusion is maximized (i.e interaction farthest from the anode) and the recoil track 
is parallel to the anode (the charge deposit is distributed along different strips). Another strong 
system requirement is to minimize the board level interconnection to allow an easy integration 
with a readout system. 



3.2 Design overview 

The front end ASIC, whose block diagram is shown in figure 0, is composed of 4 groups of 16 
channels. Each channel is composed of a current preamplifier having a gain of 15, a fast compara- 
tor (modified CMOS inverter kept in linear region) and a 5 bit DAC for setting the threshold. A 
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Figure 3. Front end ASIC block diagram. The 64 channels are decomposed in 4 groups of 16 channels. 
Each channel is composed of a current preamplifier having a gain of 15, a fast comparator (modified CMOS 
inverter kept in linear region) and a 5 bit DAC for setting the threshold. A serial link is used to configure the 
DAC. The comparator outputs are transferred serially. 



trade-off was made between the DAC resolution and the achievable preamplifier offset. The DAC 
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dynamic range, and thus its design complexity can be greatly reduced by using an autozero pream- 
plifier. This kind of amplifier measures periodically (every second) its offset during a few dozen jxs 
and determines the compensation to apply to reduce the residual output offset. The 5 bit DAC is 
designed with a 200 nA LSB (input equivalent: 13.3 nA). 

The comparator outputs are sampled at a 50 MHz rate and serialized at 400 MHz, thereby 
reducing the interconnection by a factor of 8 and as a side benefit also diminishing the power 
consumption. The serial outputs rely on the Low Voltage Differential Signaling (LVDS) standard 
to lower the electronic noise. It should be noted that using the same reference clock allows syn- 
chronous sampling between ASICs. Finally, a slow serial link is used to configure the 64 DACs 
and to individually enable/disable each channel (kill possible dead channels, ...) and to provide the 
synchronization pattern to be used. The ASIC was fabricated in austriamicrosystems BiCMOS- 
SiGe 350 nm process. It uses an effective area of 3.9 mm x 5.8 mm = 22 mm^ and requires a total 
power of 445 mW. 

4. Readout electronics 

4.1 Electronic board overview 

The readout electronic, which is an upgrade of a previous work flTDIl . comprises 8 dedicated ASICs, 
a FPGA (Field Programmable Gate Array), a flash ADC (Analog to Digital Converter) and an USB 
interface for Data Acquisition (DAQ) and slow control (see figure The connection between the 
anode located inside the chamber and the electronics at ambient pressure is done via an airtight 
interface [0]. Each strip input is equipped with a discharge protection. The hit strip information, 
that is sampled at a rate of 50 MHz, is transferred by 8 LVDS serial links at 400 MHz to a unique 
processing FPGA. This FPGA allows the auto-triggering and does the first level event building. 

In parallel to the anode signal processing, the grid signal is fed to flash ADC and sampled at 
50 MHz. While keeping a very good energy resolution, an estimation of the charge deposit through 
time can be obtained off-line by deriving the digitized CSP signal. 

The boai^d has a dimension of 25 cmx25 cm and uses 9.4 W in operation. 

4.2 FPGA firmware 

As shown in figure 0, the FPGA deserializes the data received from the ASIC and for each group 
of 16 channels a local trigger is built (OR). The first level processing starts at this stage, i.e. when 
a coincidence exists between X and Y strips, a local recording takes place (start date -i- positions). 
In theory, a single event would be defined as continuously firing strips. Unfortunately, the primary 
electron distribution can be noncontinuous, therefore untriggered strips can split the track (clusters). 
To cope with this, the recording is actually stopped when there are no more fired strips for a preset 
number of clock cycles. 

Most of the time a few strips only are fired, hence by using an adequate data encoding (shown 
in table |T]), the data payload per event can be reduced. For instance, when 2 strips in X and 2 strips 
in Y are fired in the same time slice, taking advantage of the encoding, only 64 bit are transfered 
instead of 512. This first encoding and processing stage is done in parallel for the X and Y side. At 
the following stage, dedicated state machines search and aggregate data from the same time slice in 
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Figure 4. Block diagram of the acquisition board. It comprises 8 dedicated ASICs, a FPGA, a flash ADC 
and an USB interface for DAQ and slow control interface. 
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Figure 5. Block diagram of the ASIC interface. It comprises the data deserializer, the local trigger building, 
intermediate buffering and the first level processing. 
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Table 1. Position data encoding. 



order to perform the first level event building. Tliis association is done in several stages, in order to 
concentrate more and more the data, and to present a single buffer to the USB interface. The right 
side of figure goffer a graphical representation of the FSM. Starting from the IDLE state, the FSM 
waits for the first ASIC (or group) to present data. When it is the case, the current date is saved and 
the position decoding is performed. Then the FSM looks if an other group presents position data 
marked with the same date, if yes the data are decoded and appended to the output buffer in the 
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same time slot, if not the time slot is closed and the search for a new time slot continues. At this 
stage, two possibilities remain: either no more data are available and the FSM returns to the IDLE 
state or data are still available and the earliest data has to be found. For that, a scan of the dates is 
made {seek_next_group) in order to find the first ASIC in time. This scanning is complexified by 
the fact that the time counter has a short span (13 bit which corresponds to ~ 1.31 ms), therefore it 
is performed in two stages. First the search is performed from the current date up to the maximum 
counter value, and if not successful, the current date is set to zero and the searching is performed 
again until the group is found {check roll over state). 
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Figure 6. Left figure shows the block diagram of the parallel and cascaded processing. Right figure shows 
the state machine in charge of decoding the strip coordinates and the strategy used to do the first level of 
same time slice hit coordinate data aggregation. 



The energy measurement is done in parallel to the position processing. The CSP signal record- 
ing is armed by the position triggering, but given the fact that the signal path delay is different for 
the anode and the grid signal, the actual recording is performed when the grid integrated signal is 
digitized. As a consequence of the low level signal compared to the noise level, a slope condition, 
which is more robust and noise immune than using a simple level threshold, is used. Taking ad- 
vantage of the FPGA, dedicated filtering techniques were implemented to remove low frequency 
noise (below 20 kHz) and thus to further enhance the signal to noise ratio without degrading the 
event signal shape. For that, a delayed version of the ADC signal is continuously subtracted with 
the output of a Cascaded Integrator-Comb (CIC) filter which is used to isolate the low frequency 
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noise part of the signal. Tiie resulting signal is then fed into a low pass Finite Impulse Response 
(FIR) filter in order to provide additional data smoothing. 

Consequently to the parallel processing, the position and energy data are recorded in separate 
FIFOs for USB readout. 

5. Acquisition software 

The first task of the acquisition software is to re-associate the position and the energy data. As 
shown in figure |7|, it uses the position information which is provided in a list of X/Y coordinates 
fired per time slice. Knowing that an event is defined as continuously triggering strips, the al- 
gorithm basically searches continuous triggering position (in time) and searches a discontinuity 
(time tag jumps by more than the preset value) to close the event. Once the event is defined, the 
energy information with a corresponding time tag is associated. The event building is then fin- 
ished and contains directly the coordinates of the fired strips coordinate for each time slice and the 
corresponding digitized grid signal. 

The second task of the acquisition software is to provide a real time display (see figure 
The Graphical User Interface (GUI) provides 3 projections of the track (XY-XZ-YZ), the number 
of fired strips per time slice, the Track CSP digitized signal and its derivative. It also offers the 
possibility to search (energy, duration, ...) and display a specific recoil track. An energy histogram 
of the run is also built on line and displayed on an other tab (not shown in figure ||). 

6. Summary 

A complete dedicated solution from front-end to back-end was developed to instrument a MIMAC 
prototype (256 x 256). Several MIMAC electronics can be connected per computer and the event 
building developed can be performed in several stages provided a common time tagging electronic 
board is designed and a synchronization upgrade is implemented. 
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Figure 7. Graphical representation of the event building algorithm. It uses the position information which is 
provided in a list of X/Y coordinates fired per time slice. Knowing that an event is defined as continuously 
triggering strips, the algorithm basically searches continuous triggering position (in time) and searches a 
discontinuity (time tag jumps by more than the preset value) to close the event. Once the event is defined, 
the energy information with a corresponding time tag is associated. 
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Figure 8. Screenshot of the real time display. The 3 projections of the track (XY-XZ-ZY), the number of 
fired strips per time slice, the grid digitized signal and its derivative can be seen. 



